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The magnetic excitation spectrum in the bilayer iridate Sr3lr207 has been investigated using 
high-resolution resonant inelastic x-ray scattering (RIXS) performed at the iridium L3 edge and 
theoretical techniques. A study of the systematic dependence of the RIXS spectrum on the ori¬ 
entation of the wavevector transfer, Q, with respect to the iridium-oxide bilayer has revealed that 
the magnon dispersion is comprised of two branches well separated in energy and gapped across 
the entire Brillouin zone. Our results contrast with those of an earlier study which reported the 
existence of a single dominant branch. While these earlier results were interpreted as two over¬ 
lapping modes within a spin-wave model of weakly coupled iridium-oxide planes, our results are 
more reminiscent of those expected for a system of weakly coupled dimers. In this latter approach 
the lower and higher energy modes find a natural explanation as those corresponding to transverse 
and longitudinal fluctuations, respectively. We have therefore developed a bond-operator theory 
which describes the magnetic dispersion in Sr3lr2C>7 in terms of quantum dimer excitations. In 
our model dimerisation is produced by the leading Heisenberg exchange, J c , which couples iridium 
ions in adjacent planes of the bilayer. The Hamiltonian also includes in plane exchange, J, as well 
as further neighbour couplings and relevant anisotropies. The bond-operator theory provides an 
excellent account of the dispersion of both modes, while the measured Q dependence of the RIXS 
intensities is in reasonable qualitative accord with the spin-spin correlation function calculated from 
the theory. We discuss our results in the context of the quantum criticality of bilayer dimer systems 
in the presence of anisotropic interactions derived from strong spin-orbit coupling. 


I. INTRODUCTION 

The physics of iridium-based 5 d transition-metal ox¬ 
ides has triggered considerable interest, as it represents 
the opportunity to explore the consequences of electron 
correlations in the strong spin-orbit coupling limit. Var¬ 
ious novel electronic states, topological and otherwise, 
have been predicted for the iridatesi^. The most stud¬ 
ied example to date is the relativistic Mott-like insulating 
state observed in Sr2lr04 which, in the absence of spin- 
orbit coupling, would be expected to be a metal. The 
insulating state has been argued to result from the com¬ 
bined action of strong crystal-field and spin-orbit cou- 
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pling, leading to band narrowing and a j e g = 1/2 ground 
state, on which electron correlations then act to produce 
a charge energy gap AE&. 

More generally, studies of the Ruddlesden-Popper se¬ 
ries, Sr n+ ilr„03 ra _|_i, are proving to be particularly fruit¬ 
ful, as members of this series display a striking evolution 
in their electronic and magnetic properties as the num¬ 
ber of Ir 02 -layers is increased. This can be illustrated by 
considering the end members of this series. Single layer 
Sr2lr04(n = 1 ) is an insulator (A E ~ 615 meV— ) with 
a canted, basal-plane antiferromagnetic (AF) ground 
state&i, displaying a magnon dispersion^ consistent with 
an effective Hamiltonian dominated by Heisenberg inter¬ 
actions, in agreement with theory^. Srlr03 (n = 00) 
is a strongly correlated metal with ferromagnetic cor¬ 
relations^. Bilayer Sr3lr207 (n = 2 ) represents the 
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FIG. 1. (a)-(c) Magnetic structure and excitations in 

Sr 3 lr 2 C> 7 . (a) intra and inter-layer couplings (see Eq. 1); (b) 
conventional acoustic and optical spin-wave modes previously 
used to interpret RIXS data from Sral^O?^; (c) dimer triplet 
excitation, (d)-(e) Schematic evolution of a bilayer system 
as a function of inter- to intra-layer coupling, (d) Quantum 
phase transition of the SU(2) symmetric model; (e) Effect of 
anisotropy leading to a cross over between antiferromagnetic 
order and quantum dimer regime with a gap A Es that never 
closes for finite coupling ratios. 


interesting case of being a marginal insulator (A E k, 
130 meV—), with a c-axis antiferromagnetic ground state 
(Fig.[D(a)) 12 i 13 ; and a magnetic spectrum dominated by a 
single excitation with an unusually large magnon energy 
gap A Es ~ 85 meV— . 

This excitation has been interpreted as two overlap¬ 
ping bilayer spin-wave modes (Fig. HKb) ) in the presence 
of enhanced interlayer pseudo-dipolar couplings. How¬ 
ever, in bilayer systems, the presence of a single domi¬ 
nant magnon branch is typical of weakly coupled dimers, 
in which case anisotropy generic.ally gives rise to two 
gapped modes close to each other with significantly dif¬ 
ferent intensities. Therefore, we reexamine the nature of 
the low-energy dynamics in Sr 3 lr 207 by performing high- 
resolution resonant inelastic X-ray scattering (RIXS) ex¬ 
periments which exploit a different experimental geom¬ 
etry compared to Ref. 0 . Our experiment establishes a 
fundamentally different picture of the magnon dispersion 
in Sr 3 lr 207 , with the observation of two distinct gapped 
modes. The dispersion and intensity of these modes are 
hard to reconcile with a spin-wave description, but they 
can be well accounted for by a bond-operator mean-field 
description that captures the quantum dimer nature of 
the excitations (Fig. [Qc)). 


II. EXPERIMENTAL 

RIXS is a photon in-photon out technique for the in¬ 
vestigation of the electronic structure of materials by 
probing excitations of various nature^. The scattering 
process can be described as in the following: a monochro¬ 


matic photon is resonantly absorbed by the system, pro¬ 
moting an electron from a core level into the valence 
band. This state, usually referred to as the intermedi¬ 
ate state of the RIXS process, is highly unstable and 
therefore short-lived. In RIXS one monitors the recom¬ 
bination of the core-hole to a final state of lower energy 
by a radiative transition. Analysis of the energy and mo¬ 
mentum of the emitted photon allows one to characterise 
the final state of the RIXS process; this can be either the 
ground state itself, as in elastic scattering, in which case 
the emitted photon energy coincides with the inident one, 
or an excited state. In the latter case, the energy of the 
excited state is determined by the difference between the 
incident and emitted photon energy. Beside energy, mo¬ 
mentum transfer is also used to label excitations. This 
is particularly informative when studying dispersive ex¬ 
citations, like magnetic one. 

RIXS experiments were performed on the ID20 beam¬ 
line of the European Synchrotron Radiation Facility 
(ESRF), Grenoble, with an overall energy resolution of 
25 meV. This is achieved by monochromatising the inci¬ 
dent photons with a Si(844) back-scattering channel-cut 
and using a Rowland spectrometer equipped with Si(844) 
spherical (R=2 m) diced crystal analysers^. The scat¬ 
tering plane and incident photon polarisation were both 
horizontal in the laboratory frame, i.e. 7r incident po¬ 
larization was used. The Si\ 3 lr 207 single crystal was 
grown by flux method of Ref. 112- The sample was cooled 
to a temperature of 15 K in a closed flow He-cryostat 
equipped with Be-windows. RIXS spectra were recorded 
with the incident photon energy fixed at 11.217 keV, ap¬ 
proximately 3 eV below the main absorption line as it is 
known that the intensity of magnetic excitations is max¬ 
imized at this energy^—. This shift provides a rough 
estimate of the cubic crystal-field splitting in Ir 5 d states 
consistent with previous results^. 

In Fig.[2ja), we present the in-plane momentum depen¬ 
dence of the RIXS response from Sr 3 lr 207 along high- 
symmetry directions of the Brillouin zone for L = 28.5. 
An elastic line (possibly containing contributions from 
phonons) and a magnetic excitation dominate the spectra 
in the 0—300 meV energy range. Following the conven¬ 
tion in Ref. 0, we label these features A and B, respec¬ 
tively. The latter has a sizeable dispersion of 85 ± 5 meV, 
and a gap of comparable magnitude. A weaker feature, 
C, follows the dispersion of feature B at higher ener¬ 
gies. These observations are in good accord with pre¬ 
vious RIXS measurements^. However, closer inspection 
of the dispersion along the symmetry line from (0,0) to 
(1/4,1/4), reveals the presence of an additional, previ¬ 
ously unreported feature, labelled D. This is most promi¬ 
nent around (0,0) where it is clearly separated from B and 
C. Away from (0,0), it merges almost into feature B and 
contributes to its lineshape and spectral weight. The in¬ 
tensity of feature D was found to be strongly dependent 
on the out-of-plane component of Q: it completely van¬ 
ishes when changing L from 28.5 to 25, as shown in the 
two panels of Fig. [U(b). 
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FIG. 2. In-plane wavevector dependence of the RIXS re¬ 
sponse of Sr 3 lr 207 in the 0-300 meV energy range, (a) High- 
symmetry directions, Q = (H, K, 28.5) (r.l.u.). Letters A-D 
label modes discussed in the text. White lines represent our 
model, (b) Close-up on limited region of the Brillouin zone, 
Q = (H,H, 25) (left) and Q = (H, H, 28.5) (right). 


The spectrum corresponding to Q = (0,0,28.5) is dis¬ 
played in Fig. H a )- Features B, C and D are fitted to 
three Pearson VII functions^. Feature D is clearly vis¬ 
ible, although its integrated intensity is only a fraction 
of that of B and comparable to that of C. The extracted 
dispersions of features B, C and D as a function of the 
in-plane momentum transfer for L = 28.5 are plotted in 
Fig. He). The corresponding wavevector dependences of 
the integrated intensities are shown in Fig. |4ja) . 
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FIG. 3. (a)-(b) RIXS spectra from SrsInO? for (a) Q = 

(0,0,28.5) and (b) Q = (5,0,2). (c) Energies of features 
B (red), C (green) and D (blue) as a function of in-plane 
momentum transfer. Experimental data: open circles, (H, K, 
28.5), filled diamonds, (5-H,K,2). Theory: continuous lines 
represent the transverse (red) and longitudinal (blue) modes, 
respectively. As in Fig. Ha), the dashed line is a guide to the 
eye obtained by shifting the theoretical curve for feature B by 
approximately 85 meV. 


the initial geometry. 


The results discussed so far were obtained in a geom¬ 
etry with the wavevector transfer Q predominantly per¬ 
pendicular to the IrC >2 planes. To explore how the RIXS 
spectrum depends on the orientation of Q, experiments 
were also performed with Q predominatly oriented in the 
RO 2 planes. An example of data taken in this geometry 
is shown in Fig. Hb)- It is immediately clear that the 
relative peak intensities of the features have a strong de¬ 
pendence on the orientation of Q, with feature D acquir¬ 
ing spectral weight at the expense of feature B as Q is 
rotated towards the planes. The energies and intensities 
extracted from data taken with Q in-plane are also plot¬ 
ted in Fig. He) and Fig. Ha), respectively, where they are 
seen to be in good agreement with data acquired using 


Features B and C have already been identified and 
discussed in Ref. @. Feature B was interpreted as the 
superposition of almost degenerate acoustic and opti¬ 
cal magnons (Fig. [ljb)), and their dispersion was mod¬ 
elled on the basis of a Hamiltonian which includes intra- 
and interlayer couplings, as well as pseudo-dipolar and 
Dzyaloshinsky-Moriya interaction terms. Feature C was 
assigned to the onset of a multi-magnon continuum. Fea¬ 
ture D was not reported in the previous experiment, most 
likely because of the nontrivial dependence of its inten¬ 
sity on momentum transfer both as a function of L (Fig. 
Hb)) and the orientation of Q (Fig. H a ) and (b)) re¬ 
vealed here. 
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III. THEORY 


Let us now turn to a theoretical discussion of the mag¬ 
netic excitations of Sr 3 lr 2 C> 7 . In a SU(2) symmetric 
(Heisenberg) bilayer system, there is a quantum phase 
transition as a function of the ratio of interlayer to intra¬ 
layer coupling21r— between a gapless, magnetically or¬ 
dered phase, and a gapped phase (see FigQJd)). In the 
limits of very weak or very strong inter-layer coupling, 
linear spin-wave theory and perturbation theory starting 
from isolated dimers provide very accurate descriptions 
respectively. However, these approaches fail in the in¬ 
termediate regime, and the only simple approach that 
provides a qualitatively correct description throughout is 
the bond-operator mean-field theor y 26 ’ 28- —. 

In the presence of strong anisotropy, as is the case 
in Sr 3 lr 207 , the excitation spectrum is always gapped 
(Fig. Die)), and it is impossible to know just from the 
excitation spectrum in which regime the system lies. 
However, on general grounds the excitation spectrum of 
Sr 3 lr 207 revealed by RIXS is incompatible with linear 
spin-wave theory. Indeed linear spin-wave theory pre¬ 
dicts two modes, which have dispersions related to each 
other by ca acus tic(q) — ^optical((1/2,1/2) q) (q is the in¬ 
plane momentum transfer), implying that the spectrum 
should be symmetric around (1/4,1/4) in Fig.[3])c), which 
is clearly not the case. 

We therefore developed a description of Sr 3 Ir 2 07 in 
terms of coupled dimers. In this approach, the parameter 
that controls the center of the main band is the interlayer 
coupling J c , which must then be of the order of 100 meV. 
The fact that the dispersion is approximately degenerate 
at (1/2,1/2) points to a dominant intra-plane ferromag¬ 
netic diagonal inter-dimer coupling. Finally, Hund’s rule 
exchange and the staggered rotation of the Ir-0 octahe- 
dra are known to induce anisotropic pseudo-dipolar and 
Dzyaloshinskii-Moriya interactions— We therefore con¬ 
sider the Hamiltonian 


H = J \coB(29)Su-Si J + 2wn 2 (9)Sf i Sf j + 




CiCi sin(2d)(S;i x S ij') • c z T J c ^ ) Spj • S 2 i 

i 

+ J‘2 ^ S H ■ Sij + J 3 ^ Su-Sij, (1) 


where, in agreement with Sr 2 Ir 04 , a third neighbor in¬ 
plane coupling has been included^. The naming conven¬ 
tion for the exchange constants is indicated in Fig. Oja). 
In principle, due to the staggered rotation of IrOe oc- 
tahedra, all bonds connecting opposite sublattices have 
anisotropic exchange contributions, but one can gauge 
away some of thern^, e.g. that on the interlayer cou¬ 
pling J c , which we have chosen to do. In single layer 
Sr 2 Ir 04 the angle 6 can be inferred directly from the 
canting of the in-plane ordered momenl—. In Sr 3 Ir 2 07 
the moments order along the c-axis, and 6 is a just a 



Q = (H, K, K) 



Q= (H, K, L) 


FIG. 4. Integrated intensity of features B and D as a function 
of in-plane momentum transfer, (a) Experimental data for 
features B (red) and D (blue): open circles, (H, K, 28.5), 
filled diamonds, (5-H,K,2). (b) Theory: transverse (red) and 
longitudinal (blue) responses. 


measure of the relative strength of the anisotropic inter¬ 
actions. 


A. Bond-operator mean-field method 


This model has been treated using bond-operator 
mean-field theory introduced by Sachdev and Bhatl— , 
which has proved to be an accurate method for many 
bilayer spin systems and was especially applied to spin 
ladder systems^ - — . In this formalism, first of all, a close- 
packed dimerization of the lattice is chosen. In the case of 
our bilayer system, we designate interlayer bonds on the 
sites (1, i ) — (2, *), where the first coordinate refers to the 
layer and the second one to the position inside the layer. 
The four states |tl), lit): |tt) and |44-) on each rung i 
in the Hilbert space can be combined to form one singlet 
state |si) and three triplet states |f“), a=x,y,z. To do 
so, singlet and triplet creation operators that create the 
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states out of the vacuum | 0 ) are introduced 

i*> =4io) = ^(in)-im 
i n) ^r t io) = --^dtt)-im 
\t\) ^’ t | 0) = ^= ( in ) + im 
i*?)=*f- t io) = -^(in> + iw ( 2 ) 

On each rung, the operators must fulfill a local hard-core 
constraint sjsj + Xb = 1, a = x,y, z. The action 
of the spin operators on the singlet and triplet states is 
then equivalent to 


S “i = 9 4*“ + C «/S7*f’ t *7 


js)7 


S“ — - I -«U a - t a ’U -iVf a 

°2,j — q I q s i 1 / y £g(3 r y“i q 


(3) 


i. Dispersion 

To implement the mean-field approximation, one has 
to make assumptions about the ground state. Exper¬ 
imentally, it is known that the system has easy c-axis 

I 


collinear AF magnetic order. We therefore describe the 
ground state by a condensation of singlet and triplet t z 
operators, | GS) = JX, sj, (s) ^ 0, where a new operator 
basis is defined through rotation with a rotation angle x 
adjusted to eliminate linear terms in the bosonic Hamil¬ 
tonian, 4 ( J — J 2 — J 3 ) cos(2y) = J c : 


4 = cos (x) s i ^ 
tV>t =t y ’ t, t zA = 

l l ’ % 


e*sin =t* A 

e* sin(x)st + cos(x)^’ f - 


(4) 


Here, e,; = e lC ^' Ri with q = ( 7 r, 7 r). Inserting this rep¬ 
resentation into the Hamiltonian of Eq.(l), some simpli¬ 
fications are possible: three-triplet terms have no con¬ 
tribution since they change sign under reflection along 
a plane perpendicular to the c-axis and passing through 
the centre of the rungs. Therefore, they vanish when tak¬ 
ing the expectation valued. Additionally, quartic triplet 
terms can be neglected due to their marginal effect on 
the results. This corresponds to ignoring triplet-triplet 
interactions. 

The resulting Hamiltonian can be easily diagonalized 
with the help of a Bogoliubov transformation. We ob¬ 
tain a longitudinal mode w q . z and a two-fold degenerate 
transverse mode uiq tX = u>^ y 


^q,Q: — 



with a = x, y, z and 


(5) 


Aq :2 =4J 


sin 2 ( 2 %) ( 1 — y — y 


4 J 


cos( 2 y) 


cos 2 ( 2 xbq + yb + y^q 


B„ , = 


J 


cos 2 ( 2 xbq + y< 5 q 


•h 

J 


+>q 


A q , r =2 J 


jj cos 2 (x) + sin 2 ( 2 x) f 1 “ y “ y 


\ [cos( 20 ) cos( 2 x)] 7 q + y<b + y+q, 


B ^r =w [cos( 20 ) - isin( 20 ) sin( 2 x)] 7 q + cos ( 2 xbq + cos( 2 x)<£ q , 


( 6 ) 


where <5 q = 2(cos(g x +q y ) + cos(q x - q y )), 7 q = 2 ( 0 x 357 + 
cos q y ), ifq = 2 (cos 253 , + cos 25 y ), and r = x,y. 

As long as x > 0, the gap of the transverse mode at 
q = ( 7 r, 7 r) is given by 

±e s = +y/ 4A l i +jl ±L r k 8in(29) ■ <7) 

In the absence of anisotropy (0 = 0), the model has 
two phases: i) An ordered phase with a finite staggered 
moment as long as x > 0. There is a Goldstone mode at 
(77 7 r), and this phase is thus gapless (the gap of Eq.flT]) 
vanishes for 0 = 0); ii) A gapped, disordered phase with 


no staggered magnetization when x = 0. The gap closes 
at the transition. This is illustrated in Fig. [5j(a) . In 
the presence of anisotropy ( 0 ^ 0 ), as in the bilayer iri- 
date system, there is still a phase transition at which the 
staggered magnetization disappears, but the transverse 
mode acquires a gap in the ordered phase, as emphasized 
in Eq.flTJ. Accordingly, in the disordered phase, the gap 
does not vanish at the transition. This is illustrated in 
Fig. Ob). In view of its properties, a reduced but still 
significant staggered magnetization and a large gap, we 
think that the compound S^I^CX lies in the intermedi¬ 
ate range, on the left of the transition. 
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FIG. 5. Staggered magnetization ms and magnetic gap As 
as a function of the interlayer coupling J c for (a) the isotropic 
case (9 = 0°) and (b) the anisotropic case (6 = 37°). The 
other coupling constants are equal to J = 26 meV, J 2 = 
— 15 meV and J 3 = 6 meV. In both cases, the rotation angle 
X vanishes for J c > 140 meV, bringing the system into the 
disordered phase. 



2. Intensity 


The dynamical scattering function 1^ (q, uj), /3 = 
x , y, z is proportional to the spin-spin correlation func¬ 
tion SPf( q, cej) given by 


Stf(q,w) = - ^HXlS^icfilOjfSiuj+uo-uJx) (8) 

A 


at zero temperature with the symmetric ( q c = 0) and an¬ 
tisymmetric ( q c — 7 r) rung operators defined in real space 
as Sj qo := The excited states |A) are the low¬ 

est excited states with only one triplet excitation and an 
energy uj\. The energy of the ground state |0) is ujq. Here, 
q is a two-dimensional vector describing the in-plane mo¬ 
mentum transfer and q c is the phase of the out-of-plane 
momentum of the excitations. Note that the relation be¬ 
tween q c and the out-of-plane momentum L measured 
in experiment is given by q c = 27rLd/c(mod.27r), where 
c = 20.8 A denotes the lattice parameter perpendicular 
to the bilayer, and d = 5.1 A is the intralayer distance. 

Evaluating Eq.© using the bond-operator method, we 
obtain expressions for the symmetric and antisymmetric 
part of the transverse (I q J and longitudinal (I q J inten¬ 
sities. Since the transverse dispersion branch is twofold 
degenerate, its intensity is the sum of the intensities of 
the two degenerate modes I qr = I qc + I-j c . This leads to 


the integrated intensities 

tT ( \ o-2/ \ q+7r,T + q+7i- t) 

4 (q) « 2 sin (y)-2---^-, 

q+7T,T 

jJ(q)oc2cos 2 ( X ) Aq "~ Re(J3q ’ T) , 

Wq ,r 

^(q) = o, 

1% (q) oc cos 2 (2 (9) 

Wq ,z 

Remarkably, the longitudinal intensity has no symmetric 
part. 

For a three-dimensional stacking of bilayers, the dy¬ 
namical scattering function is a linear combination of 
the symmetric and antisymmetric ones, with coefficients 
which depend on the out-of-plane momentum according 
to: 


/ T/i (q, q c ,u) 

oc cos(y) 2 /^ /L (q,w) + sin(y) 2 /J /L (q,u;). (10) 

3. Staggered magnetization 

Finally, the antiferromagnetic structure of Sr 3 Ir 2 07 
can be theoretically reproduced by calculating the stag¬ 
gered magnetization ms{ q) with q = ( 7 r, w). Since bond- 
operator theory breaks rotational symmetry in the anti- 
ferromagnetically ordered phase, the staggered magneti¬ 
zation can be obtained via the ground state expectation 
value at zero temperature of the local spin operator pro¬ 
jected along the c-axis 

ws(q) =5/Hj|(0|Sf(q)|0)| (11) 

where l = 1 or 2 refers to the layer, \xb is the Bohr mag¬ 
neton and g = 2 the electron spin g-factor. Expressing 
Eg. (fill) via bond-operator method, we obtain 

m s(q) = ^|l-a||sin(2x)| (12) 

where a = A? dk x f* dk y \v£\ 2 keeps track of the re¬ 
duction of tbe staggered magnetization due to quantum 
fluctuations. In this expression, is one of the co¬ 
efficients of the Bogoliubov transformation itq and v^, 
which can be expressed in terms of the dispersion given 
in Eq.(0) as 


Iq,z 


^q,z 


B, 




2w c 


\/^q,z(Aq, z + W q , z ) 


• (13) 


B. Discussion 

The results of the bond-operator mean-field theory ap¬ 
plied to a bilayer with significant interlayer coupling and 
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some anisotropy can be summarized as follows: i) there is 
a main band of transverse excitations whose dispersion 
reflects to a large extent the Fourier transform of the 
inter-dimer coupling; ii) the anisotropic couplings induce 
a longitudinal excitation with a well-defined dispersion 
whose energy is comparable to that of the main exci¬ 
tation band (by contrast in the SU(2)-symmetric case 
longitudinal excitations only consist of a continuum of 
two-magnon excitations); iii) the intensity of the trans¬ 
verse excitation is larger than that of the longitudinal 
one, and it peaks at (1/2,1/2), with a ratio to the small¬ 
est intensity typically in the range 1.5-3. 

Remarkably, all of the characteristics of our model 
are qualitatively consistent with the experimental results. 
Features B and D are ascribed to transverse and longi¬ 
tudinal excitations, respectively. Feature C is the lower 
boundary of the two-magnon continuum, which will be 
dominated by the dispersion of B shifted up in energy 
by the 85 meV gap, as indicated by the dashed lines in 
Figs. Ha) and[3{c). In other words, regardless of the 
details of the model, the assumption that the system 
can be described as coupled dimers with some anisotropy 
leads to predictions that are supported by our new RIXS 
data. To go beyond this qualitative observation, and in 
the absence of strong constraints provided by e.g. ab 
initio calculations, we optimized the parameter set that 
agrees best with the dispersion of both modes and their 
intensity: J = 26 meV, J2 = —15 meV, J3 = 6 meV, 
J c = 90 meV, and 9 = 37°. The dispersion curves calcu¬ 
lated for this parameter set are plotted in Fig. H a ) and 

Fig. He). 

The only aspect of the dispersion that is not accurately 
reproduced by our theory is the fact that the longitudinal 
mode seems to lie below the transverse one at (7r, 0). We 
note however that, according to a very recent improve¬ 
ment of the bond-order mean-field theory in the context 
of a 1/d expansio n 35 ' 36 , the main effect of quantum fluc¬ 
tuations for a simple bilayer model in d = 2 is to modify 
the spectrum of the longitudinal mode except at (0,0), 
while leaving the transverse mode unaffected. So this 
discrepancy is likely an artefact of the bond-order mean- 
field theory. 

The exchange pathways included in our Hamiltonian 
match those that would be obtained by projecting a Hub¬ 
bard model to fourth order—, except for the cyclic 4-spin 
terms that we have omitted for simplicity. In cuprates 
the 4-spin terms result in a zone boundary dispersion 
corresponding to an effective ferromagnetic a sim¬ 

ilar effect may explain the ferromagnetic next-nearest 
intra-layer coupling J 2 reported in S^IrO^ and also 
found here in S^I^Oy. The strong inter-layer coupling 
J c = 90 meV is qualitatively consistent with the very 
large bilayer splitting in the band structure measured by 
ARPES— Concerning the ratio of inter- to intra-layer 
coupling J c /J = 3.5, we note that in the SU(2) symmet¬ 
ric case, linear bond-operator theory overestimates the 
quantum critical ratio J c /J = 4— whereas numerical 
methods place it at 2.51—. It is therefore plausible that 



26 27 28 29 30 31 

L 

FIG. 6. Dependence of the longitudinal (blue) and transverse 
(red) intensities on the out-of-plane momentum L at the re¬ 
ciprocal space point 


treating our Hamiltonian to higher order would decrease 
the extracted ratio correspondingly by increasing J. 

With this set of parameters, adjusted to fit the disper¬ 
sion only, the intensities of the transverse and longitu¬ 
dinal modes (see Fig. Hb)) are in reasonable agreement 
with the experimental results: the response peaks around 
(1/2,1/2), and the transverse mode is more intense across 
the whole Brillouin zone than the longitudinal one. How¬ 
ever, our theory appears to overestimate the intensity of 
the transverse mode relative to the longitudinal one. 

The L dependence of the RIXS spectrum has also been 
calculated (see Fig. EO) and is qualitatively in agreement 
with the data shown in Fig. HT). Indeed, the intensi¬ 
ties vary periodically with q c , reaching their maximum 
at q c = 7r (L = 28.3) (FigH. While I T is finite for all 
momenta, I L vanishes when q c is a multiple of 27r. In 
particular, at L = 25, q c /2 tt = 6.1, i.e. q c is nearly a 
multiple of 27 t, and the longitudinal intensity nearly van¬ 
ishes, in agreement with the theoretical prediction. We 
believe that this explains why the D feature has not been 
detected in the previous experiment. 

Finally, inserting the values for the coupling constants 
and the rotation angle y = \ arccos Ja ^ j ~ 25° 

into Ea. (fl2l) . we get a staggered magnetization of 

| |sin(2x)| |1 -a\ii B = 0.76 (14) 

In this formalism, the AF order is a consequence of the 
non-zero value of the rotation angle y, which, in the 
groundstate, mixes singlets and triplets on a dimer. The 
correction due to quantum fluctuations is very small (a ~ 
0.006) because the spectrum has a large gap. As com¬ 
pared to the experimental value of about 0.52±0.08/is^, 
the theoretical value is somewhat larger. It would be easy 
to get a smaller moment by choosing a smaller rotation 
angle y at the expense of the quality of the fit of the 
dispersion. However, we have not attempted to do it 
since the bond-operator mean-field theory should not be 
considered as quantitatively accurate. 





IV. CONCLUSION 


Putting our results in perspective, we first recall that 
bilayer systems with SU(2)-symmetry have a quantum 
critical point (QCP) as a function of the interlayer cou¬ 
pling between an antiferromagnetic phase and a gapped 
phase (Fig. Old)). When anisotropy in spin space be¬ 
comes important, as in systems with strong spin-orbit 
coupling such as Sr 3 lr 207 , a new paradigm arises where 
the quantum critical point is replaced by a simple tran¬ 
sition between a gapped antiferromagnet and a gapped 
quantum dimer system, as sketched in Fig. OJe). Con¬ 
sidering the failure of linear spin-wave theory to explain 
the new mode reported in this paper together with the 
presence of a significant staggered magnetization, we are 
led to the conclusion that the system lies in the inter¬ 
mediate regime, on the left of the point where the anti¬ 


ferromagnetic order disappears, but with excitations of 
dominantly quantum dimer character. It will be reward¬ 
ing to test this conclusion with more sophisticated theo¬ 
retical approaches that could allow one to reach a fully 
quantitative agreement with experiments. 
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